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Abstract 

The Riemann hypothesis is equivalent to the Li criterion governing a se- 
quence of real constants {^k}kLii that are certain logarithmic derivatives of the 
Riemann xi function evaluated at unity. We investigate a related set of con- 
stants c n , n = 1, 2, . . ., showing in detail that the leading behaviour (1/2) Inn 
of X n /n is absent in c n . Additional results are presented, including a novel ex- 
plicit representation of c n in terms of the Stieltjes constants jj. We conjecture 
as to the large-n behaviour of c n . Should this conjecture hold, validity of the 
Riemann hypothesis would follow. 
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Introduction 

The Riemann hypothesis is equivalent to the Li criterion governing the sequence 
of real constants {\k}kLn that are certain logarithmic derivatives of the Riemann xi 
function evaluated at unity. This equivalence results from a necessary and sufficient 
condition that the logarithmic derivative of the function £[1/(1 — z)\ be analytic in 
the unit disk, where £ is the Riemann xi function. The Li equivalence [22] states 
that a necessary and sufficient condition for the nontrivial zeros of the Riemann zeta 
function to lie on the critical line Re s = 1/2 is that {\k}kLi is nonnegative for every 
integer k. 

This paper is a further contribution to our research program to characterize the 
Li (Keiper [19]) constants [22, 23]. We have previously rederived [5, 6] an arithmetic 
formula [3] for these constants, and described how it could be used to estimate them. 
Elsewhere, among several other results, we have examined summatory properties of 
the Li and Stieltjes constants, and investigated the rjj coefficients appearing in the 
logarithmic derivative of the zeta function about s = 1 [7]. In particular, a key feature 
of the sequence {7^}°^ is now known: it possesses strict sign alternation [7]. 

In this paper, we investigate a related set of constants [28] c n , n — 1,2,..., that 
might be thought of as reduced Li/Keiper constants. We show in detail that the 
leading behaviour (1/2) Inn of \ n /n is absent in c n . What remains in c n is a direct 
manifestation of fundamental properties of the zeta function. Thus, c n can be vari- 
ously interpreted as reflecting the nontrivial zeros, or the rjj constants. We present 
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additional analytic results, including an explicit representation of c n in terms of the 
Stieltjes constants 7^. We conjecture on the precise order of c n in n, we comment on 
the nature of the logarithmic derivative of the zeta function, and we briefly discuss 
possible interpretations of some of our results. 

The Gamma function is important in the theory of the Riemann zeta function- 
for instance it is needed to complete ( to the £ function. The Gamma function 
figures prominently in the functional equation for the zeta function, thereby largely 
determining the location of the trivial zeros and other analytic properties. Hence the 
digamma function appears in the logarithmic derivative of the xi function, and higher 
derivatives introduce the polygamma functions ip^) [2]. Portions of the theory of this 
family of functions are very important in much of this paper. We have also made 
extensive use of the properties of tp^ in previous works [5, 6, 7]. 

Our approach is explicit and very much in the spirit of constructivistic mathe- 
matics. Indeed, our work may be much more explicit than what would have been 
thought possible just a few years ago. 

From improved numerical calculation to height T ~ 2.38 x 10 12 [13], it is now 
known that at least the first ten trillion complex zeros of the zeta function lie on 
the critical line. For our purposes, this effectively ensures that approximately the 
first 10 26 Afc's are nonnegative, and this fact may have significant implications for our 
investigations. For instance, it may very well turn out that working asymptotically 
in k will suffice in our research program. 
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The Li equivalence is by itself a qualitative reformulation of the Riemann hy- 
pothesis. The Riemann hypothesis does not of itself dictate the exact nature of the 
Li/Keiper constants. In fact, one can easily formulate conjectures on the nature and 
order of the Li/Keiper constants that are then stronger than the Riemann hypothesis. 
These observations indicate that the Riemann hypothesis may be verifiable without 
knowing the optimal order or other properties of the Li/Keiper constants that would 
more fully characterize them. 

It is possible to use our approach also in pursuit of confirmation of the extended 
and generalized Riemann hypotheses. The corresponding A constants have been de- 
fined for Dirichlet and Hecke L-functions and other zeta functions [23], and the same 
leading behaviour 0(j In j) has been found [6]. Our attention here is strictly with the 
classical zeta function. 

Preliminary Relations 

We first recall some notation, introduce some definitions, and present an important 
Lemma for subsequent developments. We introduce the function [28] 

whose analyticity in the unit disc \z\ < 1 in the complex plane is equivalent to the 
Riemann hypothesis. Then, should the power series 

oo 

F{z) = £ c n z n (2) 

n=0 

converge for \z\ < 1, the Riemann hypothesis follows. In short order, one verifies that 
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the c n are real constants, cq = 0, and c\ = 7, the Euler constant. Therefore, we may 
write 

00 

F(z)= 1 z + Y / c n z n . (3) 

This paper investigates the behaviour of the constants c n . The importance of this 
subject is clear: from any subexponential bound on them the Riemann hypothesis 
follows [3, 28]. Indeed, in this paper we additionally conjecture the true order of 
the constants |c n |, such that this conjecture is stronger than the Riemann hypothesis 
itself. 

Suppose we knew that F(z) is analytic and univalent within the unit disc. Then 
the function F(z)/^y, satisfying F(0)/7 = and F'(0)/7 = 1, is schlicht, fulfills 
the conditions for the Bieberbach conjecture [10] to hold, and thus we would have 
\c n \ < 7n for all n = 1, 2, . . .. This shows the self consistency of the complex analysis 
involved. In fact, any direct application to the Bieberbach conjecture is thwarted due 
to the essential singularity in F, making it highly non- univalent. 

Related to a prefactor in Eq. (1), k(z) = z/(l — z) 2 is the Koebe function, and 
k a (z) = z/{\ — az) 2 with \a\ = 1 are rotations of it. These are the only functions for 
which equality holds in the conclusion of the Bierberbach conjecture. The history of 
the Bieberbach conjecture shows that it is easier to obtain results about the logarith- 
mic coefficients of a univalent function rather than for the coefficients of the function 
itself [20], and the approach of Smith [28] seems to fit within this framework. 

The classical Laurent expansion of the Riemann zeta function about the unique 
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pole at s = 1 introduces the Stieltjes constants 7& [17, 16, 25, 26], with 70 = 7. We 
have 

1 00 <—-\\ n 

C(«) = — t + E H^ 7 - (s - ir ' (4) 

ft 1 71=0 U - 

where the Stieltjes constants can be written in the form 

7 fc = lim > — In m ; . (5) 

^oo^m k + 1 J y> 

and several other forms have been given [16]. It is clear that the c„'s are multi- 
nomials in the Stieltjes constants, and that c n contains terms — (— l) n ^ n /n and 
— (— l) n 7„-i/(n — 1)!. Indeed, we are able to write much more, giving, for instance, 
an explicit formula for c n in terms of the Once again, we may therefore observe 
that sufficient estimation of the Stieltjes constants would provide verification of the 
Riemann hypothesis. 

We discuss these connections with the Stieltjes constants in a later section. For 
the moment, we simply point out that Appendix A contains explicit formulae for the 
first few c n 's in terms of them. 

The conformal map introducing the Li/Keiper constants and the coefficients c„ 
is a natural one. This map z = 1 — 1/s <-> s = 1/(1 — z) takes the right-half 
plane Re s > 1/2 to the interior of the unit circle in the complex z-plane. Just 
this sort of mapping arises in the theory of finite fields, whose zeta functions have 
zeros on a circle in the complex plane [31, 15, 4, 30]. We recall that Weil proved the 
Riemann hypothesis holds for nonsingular curves over a finite field [31], while Deligne 
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established the validity of Weil's conjectures for generalized hypersurfaces that may 
include intersections of hypersurfaces [15]. 

We first present the explicit connection between the Li/Keiper constants and the 
constants c n . We have 
Lemma 1 

— = c n + - - i hi7r + d n , n > 1, (6) 
n n 2 

where d n is the coefficient of z n of the function lnr[l/2(l — z)], and T is the Gamma 
function. That is, 

1 d n \ 1 
dn — —;-, — InT — . (7) 

We will present not only d n , but the derivatives themselves in this equation. We 
then estimate d n in n and demonstrate that the leading behaviour of \ n /n in Eq. 
(6) exactly cancels it. This is highly supportive of a decrease of \c n \ with n, and of 
the conclusion that the Riemann hypothesis should hold. Since we have previously 
conjectured as to the subdominant behaviour of the Li constants [6, 7], we thereby 
have an immediate conjecture for \c n \. 

Before proving the Lemma and going on to expressions for d n , we give some brief 
background on the Li (or Keiper) constants. The function £ is determined from ( by 
the relation [9, 11, 17, 18, 29, 27] 

^) = ^( S -i)vr- s/2 r(0c( S ), (8) 

and satisfies the functional equation £(s) = £(1 — s). The sequence {\ n }^ =l is defined 
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by 



A, 



(n-l)\ds 

The A/s are connected to sums over the nontrivial zeros of ((s) by way of [19, 22] 

l\ n 



(9) 



1-1 



and 



A 1 = - — + 1 + 1-1-2. 



(10) 



In the representation [3, 5, 6, 7] 



the sum 



— = -51 (n) + -5 2 (n) - i( 7 + lnvr + 2 ln2), 

n n n 2 



S 1 (n)= E(-l) m r (1-2— )C(m), n>2, 



rrt=2 



(12) 



(13) 



has been characterized [6]: 



n . N n 1 „ . , n . . n 1 

- Inn + ( 7 - 1)- + - < S^n) < - Inn + ( 7 + 1)- - -. 



(14) 



Further bounds on Si(n) have been developed by applying Euler-Maclaurin summa- 
tion to all orders [6]. 
For the sum 



S 2 (n) = - [ U hrn-i 



m=l 



m. 



(15) 



the constants rjj can be written as 



v i \- 1 a / k \n k+1 N s 

lim — A{m) mm 

' " k -\- 1 



fc! ^ — 1 m 



(16) 



and A is the von Mangoldt function [11, 17, 18, 29, 27], such that A(k) = Inp when k 
is a power of a prime p and A(fc) = otherwise. The constants rjj enter the expansion 
around s = 1 of the logarithmic derivative of the zeta function, 

( i^ = ~-lvAs-m k-H<3, (17) 

and the corresponding Dirichlet series valid for Re s > 1 is 

CM f-A(n) ,.„> 

The constants rjj, with ?7 = —7, have been written explicitly in terms of the Stieltjes 
constants [7, 24], a point on which we return later. Additionally, we recently proved 
the strict sign alternation of the sequence {r/j}°^ [7]. 

Now that the sum S 2 (n) has been introduced, we may present the exact relation 
Lemma 2 

^ = c. (19) 
n 

For, from Eq. (2) it follows that c n = (l/n\)(d n /dz n )F(z)\ z=0 , and from Eqs. (8) and 
(9) we have 



8=1 



-S( m )^^.M( S -DC( s )U. (20) 

Under the mapping s(z) = 1/(1 — z), derivatives transform as d/dz = s 2 d/ds, and 
the Lemma follows. 
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Proof of Lemma 1 and expressions for d n 



From Eqs. (1) and (8) we have 



ln£ 



,1 - z 
Since we have [22, 23] 



1 X = F{z) - ln(l - z) + ' 



2(z-l) 



Inn + lnT 



2(1 -z) 



ln£ 



1 



oo \ 

EA r , 



l-zj — [ n 



(21) 



(22) 



the expansion of Eq. (21) in powers of z readily yields Eq. (6). (For the latter 
equation, we have used the convention £u(z) = 2£(z), such that £l«(0) = 1, in place 
of £(z). Otherwise, Eq. (6) will have another minor term — In 2.) 

From Lemmas 1 and 2 it follows that 
Corollary 



Si 7 1 



+ 4 



n 



n 



(23) 



We next have 
Lemma 3 



dz n 



lnT 



1 



2(1 -z) 



1 d 



n-l 



2dz n ~ 1 (1-z)' 



2 (1-z) 



n+1 



l n_1 / 

1 

2(1-*) 



— l\ (n — j)! 4 



j / (1 - -2)™^' +1 d^' 



2(1-*) 
1 



2(1-*) 



^ /„ - A (n-j)! ' M (j - 1 
+ M J J(l-*)^fel 



(£-1)12^(1 -z) 



2(1 -z) 



(24) 



where ip is the digamma function and ip^ is the polygamma function. 



10 



We mention two proofs of this Lemma. The key ingredient is knowing how to 
write the successive derivatives of the digamma factor. This can be accomplished by 
applying the Faa di Bruno formula for the nth derivative of a composite function. 
We have 



dz 



2(1 -z] 



E 



/ n\ (n — 1)! 



2(1 -z] 



(25) 



(l-^) n ,t! \3J (J 1)!2'(1 ~ :H 
This equation is just a slight extension of a formula for the derivatives of a function 
9(1/ x) [8]. Equation (25), when used with the product rule, completes the Lemma. 
Another method can be based upon the expansion [1, 2] 

1 



T(z) — - exp 
z 



-~fz + 

k=2 



k z 



(26) 



giving 
lnT 



1 



2(1 -z) 



oo z n 



i n 2 k2 k 



oo 
U=0 



(27) 



Expanding the powers of power series on the right side returns us to the Lemma. 
Afterall, we have the relation 



tP U) Q = (-iy +1 j!(2 J+1 " 1)C0" + 1), J > i> 



(28) 



and this is very helpful in rewriting the constants d n below. 

From Eqs. (25) and (28) we have an exact reformulation of digamma derivatives 
of interest: 
Lemma 4 



d n 
dz n 



2(1-*) 



OO 1 
* — vn z 



z=0 



m=l 



2(1-- 

m 



n n-l 1 



2 11 



1 \"- r 

2m/ 



(29) 
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This Lemma follows by setting z = in Eq. (25), inserting Eq. (28), using the 
Dirichlet series for the zeta function, and reordering the double sum. The use of a 
derivative relation of a binomial sum completes the work. 

We are in position to write compact, yet exact, expressions for the constants d n : 
Lemma 5 



+ fit 1 



j'=i 



+ 



E 

m=l 



2 1 



m=l 
1 

m 



2 1 



m 
1 

2m 



+ 



n 



m 



2m J 
n > 1, (30) 



where ■0(l/2)/2 = —7/2 — In 2. The first line of the Lemma follows from the combi- 
nation of the results of Lemmas 2 and 3, and the second line follows from application 
of finite geometric series. 

High order approximation for the constants d n 

There are many ways in which to obtain highly accurate approximations to d n for 
large values of n. The upshot is 
Lemma 6 

1 



for 



.2(1 -z) 
That is, for j >> 1 we have 



i 1 00 

-lnvr + -^^(j)+7-ln2-l]^. 
A A .7=1 



(31) 



d; 



lnj - — - -J— + 7 -ln2 - 1 + ' 



2j 12f 



(32) 



In connection with Eq. (31), we recall the value of the digamma function at integer 
argument in terms of harmonic numbers H n : ip{n) = H n -\ — 7. 
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We indicate a couple of approaches for obtaining Lemma 6. One is based upon 
using the integral corresponding to the summation on the second line of Eq. (30). It 
turns out that this integral, h(n), was extensively studied in Appendix A of Ref. [6], 
and we have taken over the results. 

In another method, we apply Euler-Maclaurin summation to the sum over m on 
the first line of Eq. (30). In doing so, we put 



2(1-- 

m 



2 V 2m 



(33) 



such that /(l) = —2 n , /(oo) = 0, and we have the elementary integral 

1 + 2- n 



/oo 
f(k)dk = 



n 



(34) 



Therefore we obtain 



d n 
dz n 



2(1 -z) 



- 2 



-n-1 



z=0 



n 



(35) 



Then the sum over j can be performed in Eq. (30). In either approach, we discard 
any terms in the final result that are exponentially small in n, such as 2 1 ~ n /n. 

Of note, all terms in d n beyond the leading logarithmic dependence are in terms 
of integral powers of 1/n-there is no algebraic dependence upon n. 

Relations and formulae for c n 

With the aid of Cauchy's integral formula, the c n 's can be written as 

* = si/cf^ (36) 
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where C is a simple closed contour about the origin, and this can serve as the ba- 
sis of a numerical method [28]. If C is a circle of radius r about the origin, then 
c n = r~ n Jq F(re 27n<?i )e~ 27rm *<i0, and this invites the use of fast Fourier transform for 
evaluation. 

The combination of the result of Lemma 5 with Eqs. (6), (12), and (14) shows 
that c n = S2(n)/n + 0(l/n), and Lemma 2 gives the strengthening to c n = 82(11) /n. 
Since we have previously conjectured that | £2(71)1 = 0(n 1 / 2+e ) for e > 0, we have 
Conjecture 



where e > but is otherwise arbitrary. That is, we anticipate that the magnitudes 
|c n ] decrease nearly as the square root of n for large n. In Figure 1 we compare 
such a decrease with available numerical evidence [28, 24, 6]. Figure 1 contains a 
semilogarithmic plot of \c n \ versus n, together with a curve corresponding to 6/n 2 \/n. 
For this limited set, after a few initial values, the latter curve appears to provide a 
consistent upper bound. In light of the von Koch result on the Riemann hypothesis 
that tp(x) = x + 0(x 1 ^ 2 In 2 x) [14], where ip is the Chebyshev function, we suspect that 
the optimal order of \c n \ is very close to Oilnn/n 1 ^ 2 ). Figure 2 shows an example plot 
of c n versus n, that illustrates the oscillatory behaviour of these constants. In addition, 
Smith has now numerically confirmed our conjecture for the first approximately 10 5 
values of c n [28]. 

In Figure 3 we have plotted the differences 5 n = c 2 — c n _ic„ + i versus n, that 




(37) 
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appear to support a decrease in \c n \ with increasing n. In addition, the behaviour 
may indicate a correlation in the sign or other properties of the c n 's. Figure 4 plots 
the magnitude of the discrete Fourier transform applied to this sequence. This plot 
indicates underlying structure. 

Our conjecture suggests that it may be worthwhile to study in detail the properties 
of the particular polylogarithm 



oo z n 



Li/2(z) = £--17^ (38) 

n=l U 

such that L 1/2 (-l) = {^2 - l)C(l/2). 

Previously, we obtained an expression for the Li/Keiper constants explicitly in 
terms of the Stieltjes constants [7]. We recall this and related results [7]: 
Theorem 

An = l-f(lnvr+21n2- 7 )+5 1 (n)-X:(-iyf r ;) J ti £ L[% « > 2 > 

Z j=2 \J J h=l a h>0,...,j h >0 6=1 Jb- 

jiH Hh=j-h 

(39) 

5 2 (n)=n 7 -£(-iy(")j£i £ fl^, n>2, (40) 

j=2 \J / h=l U h>o,...,j h >o 6=1 

Ji+-+jh=j-h 

and 

Theorem 

Vk-i = {-l) k kj:\ £ f[ yy, k>2. (41) 

7i=l 31 >o,... Jh >0 6=1 J h - 
JiH Hh=k-h 

From Lemma 2 we obtain the exact relation 
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Theorem 



h 



7? 6 



U j=1 \3 ) h=l^ h>o,...,j h >o b=l j>>- 

JiH Hh=j-h 



n>2. 



(42) 



On the right side of Eq. (42), the constrained sum over the indices jg means that 
we have a partition of k — h over the nonnegative integers. All such partitions are 
considered, meaning that their order does not matter. The number of such partitions 
is in r?„_i or c n . 

From Eq. (4) we have 



1 -2 VI -z 



i + E 



n=0 



-In 



1-Z 



n+1 



(43) 



Then from the definition (1) and performing various expansions, we have 



Hz) = ~ E 



n=l 



n 



3=0 J 



-IV 



-li z 



\m=0 



(44) 



Carrying out the expansion in powers of z in this equation must necessarily return 
us to Eq. (42) for the coefficients c n . 

From the Hadamard product formula for the zeta function (e.g., [29]), 

exp(ln27r — 1 — 7/2)s 



C(s) 



2(s- l)r(s/2 + 1) 



(45) 



we obtain 



, ln2?r -I-7/2 , 

F(z) = ^-ln2-lnr 

1 — z 



3 - 2z 
2(1 -z) 



E ^ 



1 - 



p(l-z) 



p{\ - z) 



(46) 
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Since by Eq. (10), J2 P 1/p = °"i = ^1, we have 



F(z) 



Inn 



2(1 -z) 



In 2 -lnT 



3-2z 



1 



p{l - z) 



(47) 



2(1 -z) 

The value F(0) = I] p ln[(p — l)/p] = obtains because the sum over all the complex 
zeros of ( contains the pairs of p with 1 — p. 

From the functional equation in the form ((z) = ir z ~ 1 2 z T(l — z)((l — z) sin(7r,2/2), 
we obtain 



F(-) =F(z)-lnz + ^^ + ^— m2 + mrY— *— ] + In sin 

\ZJ Z — 1 Z — 1 VI — -2, 



7T Z 



• (48) 



.2 («-!). 

This equation should be very useful in obtaining results on the boundedness of the 



Discussion of the logarithmic derivative of ( 

This function has proved to be central in analytic number theory. Here we recall 
some known results and relate them to Eqs. (17), (18), and others. 
We have (e.g., [27, 29]) in terms of the prime counting function n(x) 



f°° ir(x) 

lnC(s) = s — -dx, Re s > 1. 

v 1 h x(x s - 1) 



x(x s — 1) 



Then 



C(s) 



s: 



tt(x) , 

-dx — s 



x s 1 ir(x) Inx , ^ 
dx, Re s > 1. 



(49) 



(50) 



C(s) h x(x s - 1) J2 (x s - If 

Since the function n has steps, these induce changes in the coefficients hence 
in ^(n) or the c n constants. In the common region of validity Re s > 1 fl \s — 1| < 3, 
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we have from Eqs. (17) and (49) 

|< S - !> + x i f ^T) & - - to < s - J ) -| ir< s - ^ < 51 > 

That is, we could write for instance 

- Hs - 1, - 1 ^ - V - K. - 1) + il f - ^ - 1) 

+ [-f^ + ^]<«- 1 >' + °K'- 1 4 (52) 

where of course by the prime number theorem ir(x) ~ rr/ In a; as x — > oo. This 
equation in powers of s — 1 gives in principle an integral representation for each of 
the coefficients rjj. 

We mention an important occurrence of the logarithmic derivative in numerical 
analysis. The reciprocal of this function is key in the classical Newton iteration for 
root finding, bringing in connections with discrete dynamical systems. Then one 
seeks the attracting fixed points of the associated Newtonian mapping. Therefore, 
from this point of view it is not unexpected that the logarithmic derivative should 
play an important role in determining zeros. 

Summary and Brief Discussion 

By way of Lemma 2, or equivalently from the use of the theory of polygamma 
functions, we have shown that the constants c n of Eqs. (2) and (3) omit the leading 
growth (1/2) Inn of X n /n. We have presented additional analytic arguments, a con- 
jecture, and partial numerical results that point to the decrease of \c n \ with n. We 
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have pointed out the limited possibility of directly applying the Bieberbach conjecture 
because the function F of Eq. (1) is not univalent within the unit disc. 

The quantity S 2 (n) is formed as the binomial sum of the alternating rjj values of 
Eq. (17). The latter is a correlated sequence. For, we have previously exhibited [7] 
the explicit summatory relation imposed upon the 77's by the functional equation of 
either the zeta or xi functions. This relation implies that a given rjj is connected to 
all the other values t/j+i, rjj+2, 

In Appendix B we call out an integral representation of the alternating zeta func- 
tion that may permit a joining of probabilisitic interpretation of the zeta function 
with Krein spectral shift functions. In turn, this may provide a useful link between 
Hardy space theory and inverse scattering theory. Though fairly independent of the 
approach of this paper, we believe it may be worth pointing this out to other inves- 
tigators. 

The importance of an explicit formula for S 2 (n) or c n should not be overlooked. 
For instance, in principle, only improved estimation of the Stieltjes constants prevents 
verification of the Riemann hypothesis by way of either the Li criterion [22] or by way 
of Criterion (c) of Ref. [3]. Concerning the magnitudes \c n \, any subexponential 
bound would serve to verify the Riemann hypothesis. 
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Figure Captions 

FIG. 1. In this semilogarithmic plot, the upper curve corresponds to values of 
6/7i 2 ^/n versus n, and the lower to values of |c n | versus n. 
FIG. 2. Plot of c n versus n. 

FIG. 3. Plot of the differences 5 n = (? n — c„_ic n+ i versus n. 

FIG. 4. Plot of the magnitude of the discrete Fourier transform of the c n sequence. 
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Appendix A: Examples of c n in terms of the Stieltjes constants 

Here, 7 is the Euler constant and 7^ are the Stieltjes constants appearing in Eq. 
(4). We have 

ci = 7» c 2 = 7-7 2 /2-7i, (-4-1) 
c 3 = 7-7 2 + ^7 3 - 27i + 77i + ^72, (A2) 
c 4 = 7 3 - \l A - ^7 2 (3 + 2 7 i) + 7(1 + 3 7 i - \l2) + I [-37i(6 + 71) + 9 7 2 - 73], (A3) 
and 

11 1 

c 5 = -7 4 + ^7 5 + 7 3 (2 + 7i) + 7p 2 (- A - 8 7i + 72) + 7[1 + 7i(6 + 7i) - 272 + -7s] 

+ ^[72 72 + 12 7l (-8 - 4 7l + 72 ) - I673 + 74]. (A4) 
The first few d^s are given by 

rf = ilnvr, (A5) 

d 1 = - 7 /2-ln2, (A6) 
d 2 = -7 + ^ 2 - 2 In 2, (A7) 

o 

4 = -3 7 + ^ 2 -61n2-^C(3), (A8) 

d 4 = -127 + ^vr 2 + - 241n2 - 21C(3), (A9) 
I lb 



and 



5 93 
d 5 = -607 + 30tt 2 + -7T 4 - 1201n2-210C(3) - -jC(5)- (A10) 
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Appendix B: The alternating zeta function in inverse spectral theory 

We first recall the alternating zeta function 

oo / i \n— 1 

I! ^ 1 — = (1 -2 1 - S )C( S ), Res>0, s^l, (B.l) 

n=l U 

this being one of the many analytic continuations of the Dirichlet series for the Rie- 
mann zeta function. Without going into the details, it turns out that this function 
can be written as an integral representation with the Krein spectral shift function 
[21] associated with the harmonic oscillator Hamiltonian on the line, with a Dirichlet 
boundary condition at the origin. As a Corollary, one may write [21] 

(1 - 2 1 ~ s )C(s) = s / e~ sx (p(x)dx, Re s > 0, (5.2) 

Jo 

where 

oo 

<l>( x ) = J2 X[ln(2n-l),ln2n](z), ( 5 - 3 ) 
n=l 

and x is the characteristic function of an interval. More generally, if the Dirichlet 
boundary condition is enforced at any other point x, a family of functions ((x, s) is 
generated. The points of discontinuity of ((x,s) satisfy a differential equation in x 
called the Dubrovin equation. This differential equation gives a curve in the space of 
analytic functions with the alternating zeta function divided by s as the initial value 
[21]. 

Now it is possible to represent the function [12] 

v (s) = ( S ~ ^ C(g) = I" e- xs Mx)dx, Re s > 0, (BA) 
s z Jo 
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with the real-valued function 



Mrf = E (1 + lnn-x). (B.5) 

l<n<e x 

Then it is possible to introduce a family of probability densities with x G [0, oo) as 
p a {x) = 4>i{x) exp(— ax)/r)(a) for a > [12]. The cumulants of p a can be written 
either in terms of the Stieltjes constants or the Li/Keiper constants at a — 1 [7, 12]. 
Comparing Eqs. (B.2) with (B.4) and (B.3) with (B.5), it appears that it should 
be possible to combine a probabilistic setting for the zeta function with an inverse 
spectral theory. This point of view offers a connection between quantum dynamics 
and stochastic processes. 
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